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^ I Abstract. Stochastic dynamics of a nonconserved scalar order parameter near its 

^^ ' critical point, subject to random stirring and mixing, is studied using the field 

^V . theoretic renormalization group. The stirring and mixing are modelled by a random 

external Gaussian noise with the correlation function ex d{t — t')k*~'^~y and the 

(^ ■ divergence-free (due to incompressibility) velocity field, governed by the stochastic 

>0 , Navier-Stokes equation with a random Gaussian force with the correlation function 

^^ ' (X S{t — t')k^~'^~y . Depending on the relations between the exponents y and y' 

jrt ' and the space dimensionality d, the model reveals several types of scaling regimes. 

Some of them are well known (model A of equilibrium critical dynamics and linear 
passive scalar field advected by a random turbulent flow), but there are three new 
nonequilibrium regimes (universality classes) associated with new nontrivial fixed 
Q I points of the renormalization group equations. The corresponding critical dimensions 

O ' are calculated in the two-loop approximation (second order of the triple expansion in 

y, y' and e — A — d). 
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1. Introduction 

Over the past three decades, increasing attention has been attracted by the dynamics 
of phase ordering — the growth of order through domain coarsening (spinodal 
decomposition), when a system (e.g. a ferromagnet or a binary alloy) is quenched 
from its high-temperature homogeneous phase into the low-temperature multi-phase 
coexistence region; see [1]-[18] and references therein. 

Much interest was focused on the late stages of the coarsening process, when some 
kind of a self-similar (scaling) regime develops with apparently universal exponents — 
the features normally associated with the critical behaviour. That regime is by now 
rather well understood; see Ref. [1] and the reviews cited the. Phenomenological 
approaches, renormalization group (RG) techniques, exactly soluble models and 
numerical simulations show that the characteristic domain size increases as a power 
of time, L{t) ~ t", where the growth exponent a depends on the global characteristics 
of the system (conserving or nonconserving dynamics, scalar or vector order parameter, 
dimensionality of space) but not on its detailed structure (like the values of the coupling 
constants). Therefore, in recent years attention has been directed to systems subjected 
to external stirring, like binary mixtures under imposed shear flow or other kinds of 
deterministic or random (e.g. turbulent) velocity flelds; see [3]-[15] and references 
therein. 

Numerical experiments and theoretical analysis (e.g. the linear stability analysis of 
the corresponding dynamic equations) of binary alloys subjected to statistically isotropic 
and homogeneous random velocity ensembles of very different kinds also suggest that, 
at least close to the critical point and under vigorous stirring, the domain growth 
is "arrested" and a new dynamical nonequilibrium steady state emerges, which is 
characterized by a continuous formation and rupture of flnite-size domains [3, 9, 16, 18]. 

Emergence of the nonequilibrium steady states appears rather a generic and robust 
phenomenon, being observed in two-dimensional numerical simulations for passive 
[9, 16, 18] and active [18] order parameters subjected to a random Gaussian velocity 
field with finite correlation length and time [9] and various kinds of regular and chaotic 
cellular fiows [16, 18]. The questions which naturally arise within this context, and which 
will be addressed in the present paper, are the following: Do those steady states reveal 
some kinds of self-similar behaviour? Do the corresponding correlation and structure 
functions exhibit power laws? If yes, do those states belong to the universality classes 
known for the models of equilibrium critical dynamics [19, 20], or do they represent 
new types of scaling behaviour? Are there any crossover dimensions for the new 
scaling regimes? Is it possible to establish the existence of these scaling regimes on the 
basis of microscopic models, and to calculate the corresponding exponents in consistent 
approximations or, better, within regular perturbation expansions? To what extent this 
behaviour is universal? What are the parameters the scaling dimensions depend on? 

We will consider the dynamics of a scalar (one-component) passive (no feedback 
on the velocity field) nonconserving order parameter ip{x) = ip{t,:K) governed by the 
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stochastic equation 

aoVtV = d'v-V'{^) + f, Vt = dt + vA, (1) 

where o"o > is the reciprocal of the kinetic coefficient and the potential V{ip) will be 
chosen as in the well-known models of critical dynamics [19, 20, 21, 22]. However, in 
contrast to the latter, the stirring noise f{x) and the velocity Vi{x) are not chosen such 
that the steady state of the system is in equilibrium, or, in other words, its equal-time 
correlation functions are not described by the Landau-Ginzburg Hamiltonian. Namely, 
the transverse (divergence-free, due to the incompressibility condition diVi = 0) velocity 
field satisfies the Navier-Stokes equation with a random driving force 

VtVi = uod\ - d{P + /„ (2) 

where V and fi are the pressure and the transverse random force per unit mass (all 
these quantities depend on x). 

The random sources f{x) and fi{x) maintain the steady state of the system and 
model the effects of external stirring and/or shaking and initial and/or boundary 
conditions. The use of such random stirring terms is a commonplace in the statistical 
theory of turbulence [23, 24, 25, 26] and other nonequilibrium phenomena [23, 27]: it 
allows one to do away with the details of the geometry of the system and to consider a 
homogeneous and isotropic problem in the infinite space. Let us specify their statistical 
properties. 

In models of equilibrium critical dynamics the form of such correlators for 
Langevin equations (like e.g. (2) without the velocity) is uniquely determined by 
the requirement that the dynamics and statics be mutually consistent (that is, the 
equal-time correlations of the dynamical problem be given by the Landau-Ginzburg 
Hamiltonian); see [19, 20, 21]. Such arguments do not apply to our non-equilibrium 
model. Like in the RG theory of turbulence, the correlators will be chosen on the basis 
of both physical and technical arguments. 

Consider for definiteness the correlation function 

D^{x, x') = {fix) fix')) = 6it - t') Dir), r = |x - x'|, (3) 

of the source field of the stochastic equation (1). The function D(r) depends only on 
r = |x — x'l, its Fourier transform being Dik). The physical arguments are that the 
noises model the injection of energy to the system owing to interaction with the large- 
scale stirring. Thus for realistic case the dominant contribution to the correlators -D(/c) 
must come from small momenta k ^ m, where m = 1/L is the reciprocal of the integral 
(external) scale L (the size of the system or a stirring device). Idealized injection by 
infinitely large modes corresponds to Dik) oc S(k). On the other hand, for the use 
of the standard RG technique it is important that the function Dik) have a power- 
law behaviour at large k. This condition is satisfied if Dik) is chosen in the form [24] 
Dik) = Do him/k) /c'^"'^"^, where Do > is an amplitude factor, d is the space dimension 
and the exponent y plays the part analogous to that played by 4 — d in the RG theory of 
critical behaviour. The function him/k) with /i(0) = 1 provides the IR regularization. 
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Its specific form is unessential; we will use the sharp cutoff h{m/k) = 6{m — k) with the 
Heaviside step function to simplify the practical calculation (in the calculations in the 
spirit of dimensional regularization one could simply set h = 1). 

The large-scale forcing is reproduced in the limit y — ^ 4, as follows from the well- 
known power-law representation of the (i-dimensional 6 function, 

6{k) = hm -i-^ / rfx {Ax)y-' exp[i(k ■ x)] = 

= S;,'k''\imU4-y){k/Ar-y], (4) 

j/^4 

where Sd is the surface area of the unit sphere in rf-dimensional space (see (41) in 
section 6) and A has the dimension of a momentum. Representation (4) also specifies 
the appropriate choice of the amplitude Dq at y — *> 4. More detailed discussion of this 
issue can be found in section 6.3 of book [21]. 

Following these ideas, the sources f{x) and fi{x) will be taken Gaussian, white in 
time (this is dictated by the principle of maximum entropy [28]) and with power-law 
spectra, (//) oc 5{t — t')k'^~'^^y for the scalar noise and (fifj) oc (5(t — t')Pjj(k)A;^"'^"^' for 
the vector one, where k is the momentum (wave vector). A; = |k| is the wave number, d 
is the space dimensionality, Pij(k) = 6ij — hikj/k"^ is the transverse projector, and y and 
y' are arbitrary parameters. The large-scale forcing corresponds to the limits y, y' — ^ 4. 

The time decorrelation of the random force guarantees that the full stochastic 
problem (1), (2) is Galilean invariant for all values of the model parameters, including 
Dq and d. As a consequence, the ordinary perturbation theory for the model (the 
expansion in the nonlinearities) is manifestly Galilean covariant: all the exact relations 
between the correlation functions imposed by the Galilean symmetry (Ward identities) 
are satisfied order by order. The renormalization procedure does not violate the Galilean 
symmetry, so that the improved perturbation expansion, obtained with the aid of RG, 
also remains covariant. 

In a wider context, the model (1), (2) can be interesting as a model system for 
studying generic nonequilibrium dynamical features. Recently, significant progress 
has been achieved in classifying large-scale, long-distance scaling behaviour of such 
phenomena, including driven diffusive systems, diffusion-limited reactions, growth, 
ageing and percolation processes, and so on; see e.g. [29] and references therein. Being 
analytically tractable, our model can serve as a good testing ground in studying such 
scaling regimes and their universality. Similar (but nonstationary) models also arise in 
stochastic inflationary models designed to describe the structure formation in cosmology; 
see [30]. 

We will apply to the stochastic problem (1) the field theoretic renormalization 
group (RG), which proved to be extremely useful in describing equilibrium critical 
phenomena, including their kinetic properties [19, 20, 21, 22]. In the RG framework, 
long-wavelength scaling regimes are associated with infrared (IR) attractive fixed points 
of the corresponding RG equations. Earlier, the RG approach was applied to the problem 
of phase separation and domain growth in a number of studies [1, 5, 6, 7, 8]. By contrast 
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with critical phenomena, however, the apphcation of the RG to the growth problems 
suffers from the lack of an (obvious) small parameter (analogous to e = A — d for 
the Landau-Ginzburg model and corresponding dynamical models), the problem also 
encountered for the stochastic Burgers and Kardar-Parisi-Zhang models [23, 27]. To 
circumvent this obstacle, in [1, 5, 6, 7, 8], the RG was used in the form of block-spin 
transformations performed using numerical Monte Carlo simulations, proposed earlier 
in [31]. Another possibility, explored in [8] (see also discussion in the review paper [1]), 
was to assume the existence of the RG symmetry and an appropriate strong-coupling 
fixed point, and then to use specific features of the conserved dynamics (absence of 
renormalization of the transport coefficient, well known for the model B of equilibrium 
dynamics [19, 20, 21]) to derive some exact relations between the critical exponents. To 
complete that analysis, however, one should take some exponents from the experiment 
or derive them using additional phenomenological considerations [1, 8]. 

The plan of the paper is as follows. We begin with the analysis of the model without 
velocity, which appears nontrivial and reveals a new type of scaling behaviour. In 
section 2 we present the field theoretic formulation of the model and its renormalization. 
After an appropriate extension, the model becomes multiplicatively renormalizable, and 
the differential RG equations can be derived in a standard fashion (section 3). The 
fixed points and their regions of IR stability are analyzed in section 4. It is shown 
that a systematic perturbation expansion in the two parameters, y and e = A — d, 
can be constructed for the coordinates of the fixed points and critical dimensions, with 
the additional assumption that y — e = 0(5^). One of the two nontrivial fixed points 
corresponds to the well known model A of equilibrium critical dynamics, while the other 
represents a new nonequilibrium universality class; the corresponding critical dimensions 
are calculated in the two-loop approximation (section 5). In section 6, the full model 
with the velocity field, governed by the stochastic Navier-Stokes equation, is studied. 
Two additional nonequilibrium scaling regimes (universality classes) are identified; the 
corresponding dimensions are found to the second order of the triple expansion in y, y' 
and e = A — d. Section 7 is reserved for a brief conclusion. Some interesting details of 
the two-loop calculation are given in Appendix A. 

The field-theoretic renormalization group was earlier applied to the problem of the 
effects of turbulence on the critical behaviour of binary mixtures in [17]. The model 
studied in that work was less realistic than the present one in two respects: turbulence 
was modelled by a Gaussian time-decorrelated statistical ensemble and the noise was 
taken to be purely thermal. On the other hand, our model is less realistic in the sense 
that the order parameter here is not conserved. Nevertheless, the main qualitative 
conclusion drawn from the two cases is the same: the instability of the equilibrium fixed 
point and the existence of a new non-equilibrium critical regime was established. Thus 
we may conclude that such a phenomenon appears quite robust and insensitive to the 
details of the model. 
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2. The model without convection: Field theoretic formulation and 
renormalization 

It is instructive to begin the analysis with the model with no convection term in (1), 
which already exhibits a nonequilibrium scaling regime and involves some interesting 
formal subtleties. The dynamical equation for the order parameter <f{x) = ip{t,:si) then 
becomes 

aodtv = 9V - V'{^) + /. (5) 

Correlator of the random noise f{x) will be taken in the form 

D^{x, x') = {fix) fix')) = 6it - t') Dir), r = |x - x'|, (6) 

with some function D(r) depending only on r = |x — x'|. The choice D(r) = 2o"o5(x — x') 
corresponds to the well-known model A of critical dynamics, which describes kinetic 
properties of the equilibrium critical state [19, 20, 21, 22]; the probability distribution 
function of its equal-time correlators is then given by exp i~Hiip)) where Hiip) = 
—cpd'^ip + Vi^p) with the implied integration over x is the Hamiltonian for the time- 
independent field v'(x). 

We assume that the model is near its critical point, and, in the spirit of the 
Landau theory, retain in Viip) only the first terms of the Taylor expansion: V"((/9) = 
TQip'^/2 + XQip^/24, where Tq is the deviation of the temperature from its critical value. 
The function D(r) in (6), however, will be chosen in the power-like form D(r) oc r"^"*"^, 
which in the momentum representation gives 

Dik) = DoA;'-'-^ (7) 

where A; = |k| is the wave number, y is an arbitrary parameter and L'o > an amplitude 
factor. The IR cutoff at A; = m is implied. 

According to the general theorem (see e.g. [21, 22]), stochastic problem (5), (6) 
is equivalent to the field theoretic model of the doubled set of fields $ = {(/?', (/?} with 
action functional 

5($) = ip'D^ip'/2 + ip' {-aodtip + 9V - roifi - Xo^Vg} , (8) 

with D^ from (6) and implied integrations over the argument x = {t, x}. Formulation (8) 
means that statistical averages of random quantities in the original stochastic problem 
can be represented as functional averages with the weight expS'(<l>). The model (8) 
corresponds to a standard Feynman diagrammatic technique with two bare propagators 
(lines in the diagrams) (v^v')o and {^p^p')o (their explicit form is given in Appendix A) 
and the vertex v^V^- 

The analysis of ultraviolet (UV) divergences is based on the analysis of canonical 
dimensions. Dynamical models of the type (8), in contrast to static models, have two 
scales, i.e., the canonical dimension of some quantity F (a field or a parameter in the 
action functional) is described by two numbers, the momentum dimension dp and the 
frequency dimension d'^. They are determined such that [F] ~ [L]^'^f[T]^'^p, where L 
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is the length scale and T is the time scale. The dimensions are found from the obvious 
normalization conditions d^ = —d^ = I, d^j^ = d'^ = 0, d^ = d'l = 0, d"^ = —d^ = 1, and 
from the requirement that each term of the action functional be dimensionless (with 
respect to the momentum and frequency dimensions separately). Then, based on dp 
and d%, one can introduce the total canonical dimension dp = dp + 2d% (in the free 
theory, dt oc (9^), which plays in the theory of renormalization of dynamical models the 
same role as the conventional (momentum) dimension does in static problems, see e.g. 
[21]. The resulting canonical dimensions are given in table 1, including the dimensions 
of the parameters which will appear later on (renormalized parameters and the others). 
It is easily checked that the role of the coupling constant (expansion parameter in 
the ordinary perturbation theory) in model (8) with correlator (6) is played by the 
combination XqDo/cxo. From table 1 it follows that this constant has the dimension A^ 
with some momentum scale A. Thus the case y < corresponds to the Gaussian IR 
behaviour (perturbation theory works in the IR range), y = is the logarithmic value, 
and for y > the RG summation is needed. The UV divergences have the form of the 
poles in y in the correlation functions of the fields $ = {(/?', (/?}. 

The total canonical dimension of an arbitrary 1-irreducible correlation function 
r = ($■■•$) is given by the relation dr = d'^ + 2d'^ = d + 2 — N^d^, where 
iV$ = {N^,N^r} are the numbers of corresponding fields entering into the function 
r, and the summation over all types of the fields is implied. The total dimension dr 
in the logarithmic theory (that is, at y = 0) is the formal index of the UV divergence. 
Superficial UV divergences, whose removal requires counterterms, can be present only 
in those functions F for which dr is a non-negative integer. Straightforward analysis 
shows that for all d > 4, superficial UV divergences can be present only in the 1- 
irreducible functions {(p'(p) with the counterterms (p'dtip, ^p'd'^ip, TQip'ip, and {(p'(p^) with 
the counterterm <f'ip^. Such terms are present in the action (8), so that the model is 
multiplicatively renormalizable. 

However, for d < 4 a new divergence appears in the function {^p'^p') (for d = 4 and 
small y, the noise correlator becomes almost polynomial in k — namely constant — 
and local in x representation). It would be erroneous to try to eliminate this divergence 
by renormalizing the nonlocal noise term (as thoroughly discussed in Ref. [32] for the 
stochastic Navier-Stokes equation). We therefore must add the local counterterm of the 
form v^V'- So we are forced to consider the variable space dimension and, for d < 4, 
to extend the original model (to include the local term in the noise correlator from the 
very beginning) in order to have multiplicative renormalizability. Finally we arrive at 
the extended model 

aodtif = (9V - T-oV^ - AoV?^/6 + /, (9) 

(fix) fix')) = 2ao6{t - t') {wok'-'-y + 1} (10) 

which has become multiplicatively renormalizable for rf < 4; the special case Wq = Q 
gives the model A (which is multiplicatively renormalizable in itself). Interpretation 
of the additional local term in (10) can be twofold. On the one hand, the fact that 
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it is generated by the renormalization procedure means that it is not forbidden by 
dimensionahty or symmetry considerations and, therefore, it is natural to include it in 
the model from the very beginning. In the language of the Wilsonian RG, this means 
that such term necessarily arises in the effective model for the properly smoothed (coarse- 
grained) field; it becomes IR relevant for d < 4, where it affects the critical behaviour and 
cannot be neglected. On the other hand, one can insist on studying the original model 
with a purely power-law correlation function. Then the extension of the model is only 
needed to ensure the multiplicative renormalizability and to derive the RG equations; 
the latter should be solved with the special initial data that correspond to the power-law 
correlator. Since the IR attractive fixed point of the RG equations is unique for any 
given choice of the parameters e and y (see section 4), the resulting IR behaviour is the 
same as for the case of the general correlation function (10) with the inclusion of the 
local term. 

By dimension the couplings Aq ~ A^ and Wq ~ A^^^ with e = 4 — d, so that we 
expect the double RG expansion in y in 5 instead of a single expansion in y (as for 
d > A) OT in e (as for the model A). However, as we will see, the real situation appears 
slightly more complicated. 

The renormalized action is 

Sr{v, ^') = a^' {wfiy~'k'~y + z^}^' 

+^' {-Z^adtV + ^3«9V - rZr,v - Z^\ii'^^/Q] , (11) 

which is equivalent to the multiplicative renormalization of the fields if -^ Z^f,ip^ 

if' — > Z' if' and parameters 

(To = Z„a, Tq = ZrT, Uq = jJ," ZuU, Wq = ^^"''Zy.W, (12) 

where we introduced the new coupling constant u = A/IGtt^ (the coefficients in RG 
functions become slightly simpler) and /i is the reference mass - additional parameter 
of the renormalized theory. 

The renormalization constants Zi = Zi{e,y,g,w) capture the divergences at 
e,y ^ 0, so that the correlation functions of the renormalized model (11) have finite 
limits for e,y = (when expressed in renormalized parameters u, w, r and /i). 

The relations between the Z's in (11) and (12) have the forms 

ZaZ^i = Zi, ZfjZ^iZ^ = 1, ZfjZ^piZ^p = Z2, 

Z^pZipi = Z3, ZT-ZipZipi = Z5, ZuZ^Z^pi = Z4. (I'JJ 

The renormalization constants Z\-Zf, are calculated according to standard rules 
from the perturbation theory; then the constants in (12) are easily found using the 
relations (13). The expansion parameter in the Z's is -u, while the dependence on the 
second coupling constant w should be calculated exactly in each order of the expansion 
in u. Like for the model A, the first nontrivial contributions to the constants Z45 
are determined by one- loop Feynman graphs, so that ^4^5 = 1 + 0(u). The leading 
contributions to the constants Z^-Zj, are determined by two- loop ( "watermelon" ) graphs 
depicted in figure 1, so that ^1,2,3 = 1 + 0{y?^. These details will be important in 
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Figure 1. Two- loop graphs giving rise to rcnormalization constants Zi-Z^. The 
simple lines correspond to the bare propagator {ff)o ~ 2a {w {k/ij,)^^^ + 1} /(w^cr^ + 
fc**) in frequency-momentum representation, whereas the lines with a slash to the 
propagator {(p(p')o = (— icrw -I- k^)~^, the slash indicating the variables of the the 
field ip'. The vertex factor is X^^. 



the analysis of the fixed points of the RG equations. The two-loop calculation of the 
renormalization constants in our model is illustrated by the example of Zi in Appendix 
A. 

3. RG equations and RG functions 

Let us recall an elementary derivation of the RG equations; see [21]. The RG equations 
are written for the renormalized correlation functions Gr = ($■■■$)/;;, which differ 
from the original (unrenormalized) ones G = ($■■■$) only by normalization and choice 
of parameters, and therefore can be equally used for analyzing the critical behaviour. 
The relation Sfi{^,e,fi) = S'(<l>, Cq) between the functional (8) and (11) results in 
the relations ^(eo, . . .) = Z^'^ Z ."^ Gr^c, fj,, . . .) between the correlation functions. 
Here, as usual, A^^ and N^i are the numbers of corresponding fields entering into F; 
Co = {o'o, tc'M^oj -^o} is the full set of bare parameters and e = {a,T,w,\ oc u} are 
their renormalized analogs; the dots stand for the other arguments (times, coordinates, 
momenta etc). We use V^ to denote the differential operation fid^ for fixed cq and 
operate on both sides of this equation with it. This gives the basic RG differential 
equation: 

{Vrg + N^l^ + N^n^,} G«(e, /i, . . .) = 0, (14) 

where "Drg is the operation V^ expressed in the renormalized variables: 

T^RG = T^,, + Pudu + (iwdw - -iaT^a - ItT^t- (15) 

In equation (15), we have written D^ = xd^ for any variable x, the RG anomalous 
dimensions 7 are defined as 

7i7 = Vfj^ln Zp for any quantity F, (16) 

and the (3 functions for the two dimensionless couplings u and w are 

Pu = i>f,u =u[-€-'ju], (17a) 

p^ = 'D^w = w[-y + e-'j,ai\, (176) 

where the last equalities come from the definitions and relations (12). 
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The anomalous dimensions 71-75 are found from the known constants Zi-Z^ (see 
Appendix A for 71). Then the relations (13) give 

7^ + 27<^/ = 7i, 7a + 27,^/ + 7^ = 0, Ja + lip' + lip = I2, 

lip + lip'=l3, It + lip + lip' = l5, lu + ^lip + lip' = l4- (18) 

Resolving these relations gives 

7a = 72 - 73, 7t = 75 - 73, lu = ll - I2 - 273 + 74, 

7^- = -71, 27^ = 73 + 72-71, 27,^' = 7i +73 -72, (19) 

and for the f3 functions (17) one obtains: 

/5„ = ^[-£-71 + 72 + 273-74], pu, = w[-y + e + ii]. (20) 

4. Fixed points and scaling regimes 

It is well known that possible scaling regimes of a renormalizable model are associated 
with the IR attractive fixed points of the corresponding RG equations. In our model, 
the coordinates -u*, w^ of the fixed points are found from the equations 

/5„(m*,w*) = 0, /3^(m*,w*) = (21) 

with the beta functions given in (20). The type of a fixed point is determined by the 
matrix Q = {Qij = d(3i/dgj}, where f3i denotes the full set of the beta functions and 
gj = {u,w} is the full set of couplings. For IR stable fixed points the matrix ^2 is 
positive, i.e., the real parts of all its eigenvalues are positive. 

From the forms of the renormalization constants (see the remark in the end of 
section 2) it follows that 74^5 = 0{u) while 71,2,3,4 = O(-u^). Therefore only 74 gives 
the leading contribution to the function j3u in (20). The actual calculation gives (see 
Appendix A) 

-fi = bu^{l+wf + 0{u^), i4 = -au{l + w) + 0{u^) (22) 

with a = 3 and b = ln(4/3). Thus the leading-order expressions for the (3 functions are 

pu = u[-e + au{l + w) + 0(m^)], 

(3^ = w[-y + e + bu\l+wf + 0{u^)]. (23) 

From (23) one immediately finds the local Gaussian fixed point u^ = w^ = 0, which 
is IR attractive for 5 < 0, e > y. The case -u* = 0, w,,, 7^ appears more subtle. 
Substituting -u* = into Pw gives f3w = w{—y + e), which suggests that for e < 0, e < y, 
the IR attractive point has w^ = 00. Then, however, the ambiguity is encountered 
in u'^{l + w)^ and the higher-order terms in j3w To resolve it, one should recall that 
for e < 0, no local counterterm to the noise correlator is in fact needed, and we can 
consider the original problem with a purely nonlocal correlator (7). For such a model, 
one easily finds that y < corresponds to the Gaussian fixed point with finite amplitude 
of the noise correlator (in fact, for the purely nonlocal case in can simply be scaled out 
from the action). The case y > 0, e < requires a bit more elaborate analysis, which 
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shows that the large-scale behaviour also in this region is that of the Gaussian model 
with the nonlocal correlator. Thus, we finally may conclude that the region e < 0, 
e < y corresponds to the Gaussian fixed point with -u* = and purely nonlocal noise 
correlator. 

In models with two regulators like e and y, it is usually implied that they are of 
the same order, y = 0{e), and the coordinates of the fixed points and the values of 
anomalous dimensions at those points are sought in the form of double series in e and 
y; see e.g. [32]. No such solution, however, can be constructed in the case at hand, due 
to the absence of an 0{u) term in the square brackets for f3^ in (23). A regular solution 
can be constructed if we assume that y = 0{e), but their difference y — e is of order 
e"^. In other words, the actual expansion parameters appear to be e and ^y{e — y). It 
is reminiscent of the well known y/e expansion for the Ising ferromagnet with quenched 
disorder [33] . There it was a consequence of an accidental degeneracy of the P functions 
of the two involved couplings, which for their ratio (the analog of our w) also implies 
vanishing of the first-order contribution. A similar situation was also encountered in 
[17] where the mixing of a conserved order parameter by a Gaussian velocity ensemble 
was studied. 

To be definite, let us write y = e + Be^ with some B. Now we have 

13^ = u[-e + au{l + w) + 0(m^)], 

p^ = w[-Be^ + hu^{l + wf + 0{u'^)] (24) 

and the solutions for m* and w* can be found as regular series in e, while the dependence 
on B should be taken into account exactly in each order of the e expansion. Now we 
can identify two nontrivial fixed points, which we denote as I and II. 

For the point I, we find w^, = 0, -u* = e/a. This point clearly corresponds to the 
model A of critical dynamics; w* vanishes identically to all orders of the e expansion (the 
model with wq = w = Q is local and therefore closed with respect to the renormalization), 
while M* has nontrivial corrections of order e^ and higher. The Vt matrix at this point 
is: 

Puu 
Puw 
Pwu 

(3^^ = -Bs^ + hu^ = -Be^ + bie/af. (25) 

It is triangular, so the fixed point is IR-attractive if the diagonal elements f3uu and f3ujw 
are positive. This gives e > (of course) and B < b/a^ . Thus we have established that 
the region of IR-stability of the fixed point I includes not only the sector 5 > 0, e > y, 
but is slightly wider: it also involves a narrow "beak" adjacent to the ray e = y in the 
region e < y; see figure 2, where regions of stability of the nontrivial fixed points of the 
model (9), (10) are depicted. 

Furthermore, we can see that the fixed point with both u^, ^ and w* 7^ also 
exists for the beta functions (24); we denote it as point II. Indeed, Pu, = and -u* 7^ 



a-u*(l 


+ w,) = 


'■£, 




aul = 


= eVa, 






2u^bw^.{l + w 


*)- 


= 0, 


-Be 


:' + bu'-- 


= - 


-Be' 
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Figure 2. Regions of stability of the fixed points of the extended model (9), (10). The 
region labelled GNL (Gaussian nonlocal) corresponds to the fixed point with vanishing 
interaction {d > 4) and purely nonlocal correlation function of the noise; the label NL, 
u 7^ refers to the nonlocal regime II of the extended model; the first-order boundary 
between the latter and the region of stability of the universality class of model A 
(the local regime I of the extended model) is y — e and the improved boundary is 
y ^ e + BcS^ . The boundary between the universality class of model A and that of 
the linear model driven by the random noise with local correlations only - labelled GL 
(Gaussian local) - is the negative y axis. 



gives u^{l + w^) = e/a, substituting in /9^^ = and assuming w^ j^ gives 

= -Be^ + 6(1 + w,)K(l + w,)f = e^-B + b{l + w,)/a^} , 

so that at the fixed point 

w* = -1 + Ba'^/b, M* = eb/Ba^ 

with corrections of order 0{e) and 0{e'^) respectively. 
The Q matrix at this point is: 

f3uu = au^{l + w*) = e, 

Puw = aul = £^6V^^«^ 

/5^„ = 2u^bw^{l + w^f = 2bw^{e/a) {a^B/bf, 

(3ww = Sw^bul{l + w^y = 3bw^{e/a)'^. 

Thus the Q matrix has the form 

a As 



(26) 



n = e 



(3 Be 



(27) 
(28) 
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Although the elements in the right column are small in e in comparison to the left 
column, they are needed to find the eigenvalues to leading order. It is important that 
the 0(6^) corrections to the left elements are not needed (they give only corrections), 
despite the fact that they would be of the same order as the right elements. 
The eigenvalues are: 

Q^ = a6 + 0{e^), Q2 = {B-f3A/a)e^ + 0{e^). (29) 

The point is IR-attractive if f2i 2 > 0. Substituting (27), (26) gives: e > 0, B > h/a^. 

Thus the curve y = e + BcS^ with B^ = h/a^ is the boundary between the regions 
of IR-stability for the point I (e > 0, -B < h/a^) and the new nonlocal point II {e > 0, 
B > h/a^). There is neither gap nor overlap (at least in this approximation). 

The physics requires that the coordinates of any physical fixed point, -u* and w*, 
be non- negative {u and uw are amplitudes in pair correlators). One can check that 
this condition is automatically satisfied in the region of their IR stability. For example, 
w* > in (26) gives B > b/a"^, which is also ^^2 > 0, and so on. 

The resulting pattern of the regions of stability of the nontrivial fixed points of the 
model (9), (10), shown in figure 2, is as follows: the quadrant e > 0, y > is divided 
into two parts by the parabola y = e + B^e'^ with B^. = b/a^. The part below it is the 
region of IR stability of the point I (universality class of the equilibrium model A); it 
includes the whole sector e > Q, e > y. The part above it is the region of IR stability 
of the point II. It corresponds to a new nonequilibrium universality class, where the 
nonlocal term in the random force is important. It is worth noting that Bf. appears 
rather small: b = ln(4/3) ^ 0.287683 and B^ = ln(4/3)/9 ^ 0.032. 

5. Scaling behaviour in the IR range 

Existence of IR-attractive fixed points implies scaling behaviour with definite critical 
dimensions A^ of all quantities F (fields and parameters): 

Ai. = 4 + A^rf^ + 7^, A^ = 2 + 7:, (30) 

where dp'^ are the canonical dimensions of F, given in table 1, and 7^ is the value of 
Jf at the fixed point in question; see e.g. [21, 26]. In the case at hand 7^. = 'yp{u^,w^). 
In particular, for the pair correlator of the field ip this gives: 

(V9(x + r, t + t')v^(x, t')) = r-^^- ^ (tq r^^ t/r^-) . (31) 

This is the leading term of the asymptotic behaviour in the IR range, determined by the 
inequalities /c ~ 1/r ^ A, where A is the UV momentum scale defined by the relation 
■uq ~ A^. JF is a universal scaling function of two arguments, which are supposed to be 
of order unity; this completes the definition of the IR range: tq r^^ ~ 1, tjr^" ~ 1. (In 
the free theory we would have ror^ ~ 1, t/r^ ~ 1). The correct canonical dimensions 
in (31) are guaranteed by the amplitudes built from IR irrelevant parameters a and A, 
not shown explicitly. One usually assumes that T has finite limits for tq = (that is. 
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Table 1. Canonical dimensions of the fields and parameters; d = 4 — e. 

F ip if' Vi v'i (To, a uo,u m,/i,A r Aq wq go \,w,g 

4f-lf-l-lrf + 12 -2 1 2 4:-dd-4: + yy'0 

d"^ 1 1-1-110 

dirf-lf + 11 d-1 1 2A-dd-A + yy'0 



exactly at the critical point) and/or t = (equal-time correlator). Then from (31) one 
obtains 

(V9(x + r,t)v9(x,t))=r-2^-. (32) 

From the general expression (30), canonical dimensions dp^ from the table 1 and 
the relations (19) we obtain 

A^ = 1 - e/2 + (73 + 72 - 7i)/2, A^ = 2 + 72 - 73, 
A^, = 3-5/2 + (71 + 72 -73) A A, = 2-73 + 75- (33) 

In the leading approximation the anomalous dimensions have the forms 

-fi = bu^{l + wf, 72 = 6M^(l + w)^ 73 = m2(1 + w)V6, 

74 = -3m(1 + w), 75 = -m(1+w) (34) 

with b = ln(4/3). For w = these expressions coincide with the results known for the 
model A; see e.g. [21]. 

Consider first the fixed point I with -u* = e/a + O(e^) and w^ = 0. The local 
model {w = 0) does not generate nonlocal counterterms and is "closed with respect 
to renormalization." As a result, w* vanishes identically, the dimensions (33) are 
independent of y and coincide with their analogs of the model A to all orders of the e 
expansion. For the model A, the fluctuation-dissipation theorem gives the exact relation 
7^ = 72, so that these quantities disappear from the expression for A^^ in (33), while 
the dimensions 'yt.ru (and hence A^^^) coincide with their counterparts in the static Xip^ 
model; see e.g. [21, 22]. 

The standard notation for the equilibrium case is 

A^ = d/2-l + r]/2, A, = l/z/, A^ = z, (35) 

while for (/?' the above identities give A^/ = A^ + z = d/2 — 1 + z + ri/2. The 
static exponents rj and u are well known from 4 — e, 2 + e and 1/A^-expansions, 
real-space RG (all augmented by various summations), high-temperature expansions 
for the Ising model (considered most reliable), Monte-Carlo simulations. The values 
recommended by [21, 22] are r/ = 0.0375 ± 0.0025 and v = 0.6310 ± 0.0015 (Borel 
summation of 5-order results). For z only two terms of the 4 — e expansion are known: 
z = 2 + 0.726(1 — 0.1885e)?7 [34]; there are also leading-order results in 2 + e and 
1/A^-expansions (see the references in [19, 20, 21]). 
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To avoid possible misunderstandings, it is worth noting that the scahng behaviour 
of the model A and the extended model (11) near its fixed point I coincide only to the 
leading orders given by the expressions (31), (32). The corrections to those expressions 
are different. In particular, the leading correction from the UV range has the form 
[l + Ci{k/A)^^ + C2{k/A)^'^~\, where f2i^2 are the eigenvalues of matrix f2 for the fixed 
point I given in (25). One of them involves the parameter y which is absent in the pure 
model A. 

Let us turn to the nonlocal regime of the extended model (11), described by the 
fixed point II. Substituting (34) and (26) into (33) gives 

A^ = 1 - e/2 + (£V2)(l/54 + b/9 - B) + 0{e''), 

A^, = 3 - e/2 + (£V2)(-l/54 + b/9 - B) + 0{e^), 

A^ =2 + (5/3)2(6- 1/6) + 0(5^), 

A, =2-e/3 + 0{e^). (36) 

We recall that B comes from the relation y = e + Be^ and the nonlocal fixed 
point is IR-attractive when B > b/9 ~ 0.032, see (29). Note that the dimension 
A^ = 2 + 0.00136:2 + 0{e^) appears surprisingly close to its canonical value A^^ = 2. 

Consider the real case d = 3, then y = 1 + B and the noise correlator is l/r^~^. 
The exponent 77, defined by the same "equilibrium" relation A^p = d/2 — 1 + ri/2, takes 
on the form rj = (1/54 + b/9 — B) ^ 0.05 — B. So rj can be made negative for reasonable 
B: in particular, for i? = 1 and noise correlator l/r^ we have 77 ^ —1. In this respect, 
the nonequilibrium steady-state scaling differs from the equilibrium case, described by a 
local Landau-Ginzburg action: for the latter, the exact inequality rj > can be derived 
from the unitarity of the corresponding pseudo- Euclidean quantum field theory [35]. 
Bearing in mind possible cosmological application of the model [30], it is tempting to 
note that rj = —1 corresponds to the Zeldovich spectrum [36]. 

6. Inclusion of the velocity field 

Let us turn to the full stochastic problem (2), (9), (10). The field theoretic action 
functional then becomes 

Si<^) = S,{v',v) + if'D^ip'/2 + ip'{-aoVtV + ^V - ^oV? - W/6} , (37) 

where 

S,{v',v) = v'DW/2 + v' {-Vt + i^od^} V (38) 

is the action functional for the stochastic problem (2), D^p and D^ are the correlation 
functions of the random forces / and /,, respectively, Vj = dt + vidi, and all the required 
integrations over x = {t, x} and summations over the vector indices are understood. The 
new full set of fields $ = {ip,ip',Vi,v'j} involves the auxihary vector field v'^. It is also 
transverse, div'^ = 0, which allows one to omit the pressure term on the right-hand side 
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of relation (38). Correlation function D^ is given by (10), while D^ will be taken in the 
form 

{Mx)f,ix')) = govl 5{t - f)k''"'~y'R,,ik), (39) 

where -Pij(k) = 6ij — kikj/k"^ is the transverse projector, y' a new arbitrary parameter 
analogous to y from (10), go ~ A^ is a new positive coupling constant, and the factor z/q 
is explicitly isolated for convenience. The IR cutoff ai k = m is also implied. Canonical 
dimensions of all the new parameters and their future renormalized counterparts are 
given in table 1. 

The stochastic Navier-Stokes equation (2) with a power-law noise spectrum was 
introduced a long ago [23, 24, 25] and is by now very well studied, at least for small 
values oiy'. The two-loop results have been derived recently in [37]. Detailed exposition 
of the RG approach can be found in [21, 26]; below we confine ourselves to only the 
necessary information. 

The model (38) is logarithmic (the coupling constant go is dimensionless) at y' = 0, 
and the UV divergences have the form of the poles in y' in the correlation functions of 
the fields v, v'. Dimensional analysis, augmented by some additional considerations 
(Galilean symmetry and structure of the vertex), shows that for all d > 2, the 
superficial UV divergences, whose removal requires counterterms, are present only in 
the 1-irreducible function (v'v) , and the corresponding counterterm reduces to the form 
v'd'^v. Owing to the form of the vertex (the derivative can always be moved onto v' 
using integration by parts), the divergence in the function {v'v') (allowed by dimension 
for (i < 4) is in fact absent for all d > 2. So the local counterterm v'v', analogous to (/?'(/?' 
in (11), is not needed here. For this reason, we did not include the constant contribution 
to (39), in contrast to its scalar counterpart (10). Then for the complete elimination of 
the UV divergences it is sufficient to perform the multiplicative renormalization of the 
parameters z/q and g^ with the only independent renormalization constant Z,^: 

VQ = VZy, go = gfJ^'^'Zg, Zg = Z^^. (40) 

Here /i is the reference mass in the MS scheme, g and u are renormalized analogs of the 
bare parameters go and z/q, and Z = Z{g,y',d) are the renormalization constants. In 
contrast to the model (11), no renormalization of the fields is needed, Z^,^r = 1. The 
relation between the Z's in (40) results from the absence of renormalization of the noise 
term in (38). Now the standard RG equations are readily derived, the corresponding 
function f3g = V^g in the one-loop approximation is 

PM = 9{-y' + 37.) = 9i-y'+ 5f~wf'f. l + 0{g') (4i) 



4(rf + 2)(2 



TT 



where Sd = 2%'^^'^ /T(d/2) is the area of the unit sphere in d dimensions. From (41) we 
immediately conclude that for y' > 0, the model has an IR attractive nontrivial fixed 
point (7* > iPgig*) = 0, P'g{g*) > 0), while for y' < the IR attractive fixed point 
is Gaussian, g^, = (that is, the nonlinearity in (2) is IR irrelevant). From the first 
equality in (41), which follows from the relation between the Z's in (40), the value of 
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'-fi, at the nontrivial fixed point is found exactly: 7* = y'/3 (no corrections of order {y')^ 
and higher). As a result, critical dimensions of the frequency and the fields are also 
found exactly: 

A^ = 2-y'/3, A, = l-y'/3, A,, = d - I + y' /3. (42) 

These results remain intact in the full model (37), because the inclusion of the 
scalar fields does not affect the velocity (the field ip is "passive"). We thus conclude 
that for y' < 0, where the IR behaviour of the velocity becomes Gaussian, the scaling 
regimes of the full problem (37) are described by the fixed points I and II from section 4 
depending on the relation between y and e = 4 — d. 

For y' > the contributions of the velocity field become IR relevant, and the RG 
analysis of the full problem (37), (38) is needed. For d > A, where no local counterterm 
ip'ip' is required, the model in question is formally equivalent to the stochastic model of 
the turbulent advection of a chemically active scalar field, studied earlier in [38]. We 
will be interested in the case d < 4, where the local counterterm should be included from 
the very beginning to ensure multiplicative renormalizability. Then the UV divergences 
manifest themselves as poles in the full set of regulators y, y' and e = 4 — d. Dimensional 
analysis and symmetry considerations show that, for all d > 2, the full model is 
multiplicatively renormalizable, and the corresponding renormalized action has the form 

+ V?' {-Z2aVtV + Zg^V - rZ,^ - Z^Xfi'^ye] , (43) 

where Syji{v',v) is the renormalized analog of the action (38) expressed in renormalized 
variables using relation (40). The renormalization constants Z\-Zf, now contain the 
poles in the full set of regulators y, y' , e = 4 — d and, in comparison to (11), depend 
on two additional dimensionless couplings g and s = au. One can easily see that 
the relations (12), (13) for the new constants Zp and (18), (19) for the corresponding 
anomalous dimensions '-fp remain valid in the extended renormalized model (43). The 
RG equation takes on the form 

{Vrg + N^l^ + N^n^, + Ar,7, + Ar,,7,,} ^^(e, /i, . . .) = 0, (44) 

where I^_rg is the operation D^ expressed in the renormalized variables: 

Vrg = V^ + f3gdg + pudu + Mu, + Psds - l.V, - ^^V,, (45) 

with Pg from (41), Pu,w from (17), and the new function Ps is 

Ps = -S7s = -sIJu + la] = -sl'ju + 72 - Ts] • (46) 

Furthermore, it is easily checked that in the leading order, the anomalous dimensions 
7i, 2,4,5 remain the same as in the model without velocity and can be taken from (34). 
On the contrary, 73 acquires an additional 0{g) term, in comparison to which the O('U^) 
term in (34) is only a correction and should be neglected. So the leading-order expression 
for 73 becomes 
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As usual, scaling regimes of the full model are associated with the IR attractive fixed 
points, whose coordinates are found from the equations jS^ = with x = {g,s,u,w}. 
Due to passivity of the field 93, the function Pg is independent ofx = {s,u,w} and 
the corresponding elements df3g/dx of the matrix Vl vanish. Thus Vt is triangular, its 
elements dPx/dg with x = {s, u, w} do not affect the eigenvalues, and we can set g = g* 
in the functions /5^ with x ^ g from the very beginning. Now we are treating e = A — d 
as one of the small expansion parameters, and in the leading-order approximation we 
should set (i = 4 in (41) and g^. In the leading approximation, the function I3s is 
independent of u and w, and the elements dfis/dx with x = {u,w} also vanish. The 
element d(3s/ds coincides with one of the eigenvalues of Q. It is easy to see that, for 
y' > 0, this eigenvalue d^s/ds cannot be positive for s^, = 0, so that the condition /^^ = 
implies 7;/ + 72 — 73 = 0, which in our approximation gives 3s^ — s* — 1 = 0. Physical 
considerations require s* > 0, which finally gives s,, = (1 + vT3)/6 ~ 0.7676. Since the 
elements df3x/ds with x = {u,w} are again irrelevant in the analysis of the remaining 
eigenvalues, we can set s = s* in Pu,w We are left with the system of two equations, 
Pu,w = 0, where 

f3u = u[-e + au{l + w) + 0{g'^,u^)], 

(3^ = w[-y + e + bu\l + wf + 0{g\ u% (48) 

Here we have denoted 

e = e-2y'/?,, y = y-2y'/?,. (49) 

Up to the notation, this system coincides with that describing the model without 
velocity, equation (23), and its fixed points and their regions of stability can immediately 
be inferred from the results of section 4. In order to have a regular expansion we again 
set y = e + Be'^, while e and y' are of the same order. (From the relation e — y = e — y one 
obtains Be^ = Be^. We could also write y' = B'e, but this is not necessary.) Now we 
can identify two fixed points, which we denote III and IV, corresponding to completely 
new universality classes. 

The point III corresponds to a purely local model: w^, = to all orders of the 
perturbative expansion, and -u* = e/a with corrections 0{e'^). It is IR attractive for 
e > and B < h/a^ (we recall that a = 3 and h = ln(4/3)). For the point IV, one 
obtains w* = — 1 + Ba'^/b and m* = eb/Ba^, with corrections of order 0{e) and 0(£^) 
respectively. It is IR attractive for e > 0, B > b/a^. The curve y = e + B^e^ with 
Be = b/a^ ~ 0.032 is the boundary between the regions of IR-stability for these two 
points. 

For both these regimes, the relations /3t, = and f * 7^ imply 

i: = ii-ii = -i: = -y'i^. (50) 

and the critical dimension of frequency is found exactly: 

A^ = 2 + 7: = 2-7: = 2-|/V3, (51) 
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in agreement with (42). Combining Eqs. (19) and (50) gives 

27; = 73* + 1*2 - 7i = 7: + 272 - 7i*, 

27;' = 73+7i-72* = 7: + 7i, (52) 

which, along with (51), allows one to determine the dimensions A^ = {d/2 — 1) + 7* 
and A^r = (^d/2 — 1) + A^ + 7*, in the second order of the e expansion without practical 
calculation of the two- loop corrections to 73 and v^: 

A^ = 1 - e/2 + y'/6 + be^/18, A^/ = 3 - e/2 - y'/6 + be^lS (53) 

for the fixed point III and 

A.^ = 1 - e/2 + y'/6 + be^/9 - Be'^/2, A^, = 3 - e/2 - y'/6 + Be^/2{5A) 

for the fixed point IV, with cubic-in-e corrections. (We recall that in counting the orders 
we imply e ~ e ~ y', e — y ~ e^). In fact, the expression for A^^r in (54) holds to all 
orders (no corrections of order e^ and higher), because '-fl = y — e = Be^ exactly, as a 
consequence of the relations /9^ = and w^ 7^ 0. Finally, for the both points III and IV 
one obtains A,- = 2 + e/3 + y'/9, with quadratic corrections. 

For e < 0, the self-interaction of the scalar field becomes irrelevant (-u* = 0), and 
we obtain two more fixed points, which correspond to the scalar field subject to a linear 
diffusion-advection equation, with the velocity ensemble given by the action Sy{v',v) 
from (38). The regions s > y and s < y correspond to the local (w* = 0) and nonlocal 
correlator of the scalar noise, respectively. For the purely nonlocal case, the amplitude 
w remains finite at the fixed point and can be eliminated from the action by appropriate 
rescaling of the fields ip and (/?'. 

The linear passive scalar case is very well understood within the RG framework 
(see e.g. chapter 2 in the book [26] and references therein). The only superficially 
divergent 1-irreducible correlation function is {(f'(f) with the counterterm ip'd'^^p, and 
the corresponding renormalization constant (in our notation identified with Z^p = Z^^) is 
independent of the form of the noise correlator (the latter only determines the canonical 
dimensions). The absence of renormalization of the noise term results in the exact 
relation Z^/ = 1 (it is implied that the amplitude is scaled out from the action). Along 
with the relations (50) and (51), which remain valid for the passive linear case, this 
gives exact results for the dimensions: 

A^ = d/2-l + y'/3, Ap, = d/2 + l-y'/3, A, = 2, (55) 

which in the standard "equilibrium" notation corresponds to rj = y'/3 and u = 1/2, 
different from their counterparts for the standard model A. 

7. Conclusion 

We have studied stochastic model that describes dynamics of a nonconserved scalar field 
(order parameter) near its critical point, subject to random external stirring and mixing, 
in d spatial dimensions. The stirring was modelled by an additive random Gaussian noise 
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with the pair correlation function oc 6{t — t')k^~'^~y. The mixing was modelled by the 
convection term with a divergence-free (due to the incompressibility condition) velocity 
field, governed by the stochastic Navier-Stokes equation with a random Gaussian force 
with pair correlation function oc S{t—t')k^~'^~y . Possible scaling regimes of the model are 
associated with nontrivial IR attractive fixed points of the corresponding RG equations. 
Their coordinates, regions of stability, and the corresponding critical dimensions can be 
calculated within a systematic expansion in y, y' and e = 4 — rf (or only y and e = A — d 
for the model without velocity) with the additional assumption that y — e = 0{e'^). 
Depending on the relations between those parameters, the model reveals several types 
of scaling regimes. Some of them are well known: model A of equilibrium critical 
dynamics and linear passive scalar field advected by a random turbulent flow, but there 
are three new nonequilibrium universality classes, associated with new nontrivial fixed 
points. In this sense, the critical behaviour of the model appears richer and less universal 
than that of the equilibrium critical dynamics. 

The critical exponents (dimensions) for the new universality classes are derived in 
the second order of the expansion in y, y' and e (two-loop approximation). 

It remains to note that the large-scale mixing {y = y' = 4) in three dimensions 
(6 = 1) belongs to the universality class of the linear passive scalar with the nonlocal 
noise correlator and therefore corresponds to the dimensions (55). Of course, the results 
of our perturbative RG analysis are absolutely reliable and internally consistent only 
for small values of the expansion parameters e, y and y', while the possibility of their 
naive extrapolation to finite (and not small) real values is far from obvious. On the 
other hand, the observation that the (y^'^-interaction becomes irrelevant for the large- 
scale forcing is reminiscent of the results derived in [13, 14]. There, it was argued that 
a non-random shear flow strongly suppresses critical fluctuations, and the behaviour of 
the system becomes close to mean field in the strong shear limit; see also discussion in 
[15]. 

Our analysis can be directly generalized to the cases of a A^-component order 
parameter, presence of anisotropy, compressibility etc. The generalizations are 
straightforward but rather cumbersome (for the stochastic Navier-Stokes equation, see 
e.g. chapter 3 in the book [26] and references therein). On the contrary, the case of 
a conserved order parameter appears rather different from both the conceptual and 
technical viewpoints (namely, it involves two different dispersion laws: cj ~ /c^ for the 
velocity and a; ~ /c*^ for the scalar). These issues will be addressed elsewhere. 
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Appendix A. Calculation of the renormalization constants 

Consider as an example the calculation of the constant Zi and the anomalous dimension 
7i for the model (10) without the velocity field in detail. The leading contribution here 
is given by a two-loop Feynman graph, so this example is representative: calculation of 
the other renormalization constants (including the velocity field) can be performed in a 
similar way (for two-loop contributions) or is much easier (for one- loop graphs). 

The 1-irreducible function P = (v^V') i^ the renormalized critical (r = 0) theory to 
order O(A^) has the form 

\2 2e 

T = 2a \w(k/fxy-' + ZA + —^— x Diagram, (A.l) 

I. J g 

where the first term is the noise correlator written in renormalized variables, Zi should 
be taken to order O(A^), 1/6 is the symmetry coefficient, A^yU^^ comes from the vertex 
factors, dashed external ends correspond to fields (/?' and the three identical lines 
correspond to bare propagators 

2a {w jk/fir-^ + 1} 

W'V'/O = 2 2 1/4 '^^•2) 

in frequency-momentum representation and 

{^ip)o = {w ik/fi)y-' + 1} 1 exp (-P|t - t'\/a) (A.3) 

in time-momentum representation. The second bare propagator 

(y9V9')o = {-i(^uj + A;2)-i ^ e{t - t') exp {-k\t - t')/a) (A.4) 

existing in the diagrammatic technique for the model (11) does not appear in this 
diagram. Within our accuracy, all Z's occurring in the diagram (for example, Z| coming 
from renormalized vertices) have been replaced with unities. 

The constant Zi will be found from the condition that it cancels the poles in e and 
y which are present in the diagram, so the full expression (A.l) is regular in e and y and 
finite for 5 = 0, y = 0. This requirement determines Zi up to a regular part. We will 
use the minimal subtraction (MS) scheme where Zi = 1+ only poles in e, y and their 
linear combinations. In the model with two regulators like e, y there are subtleties in 
defining the MS scheme in higher orders (for example, is e/y a pole or a finite thing). 
We will work only in the leading order and can neglect these subtleties. 

It is sufficient to calculate the diagram with the external frequency and momentum 
equal to zero. In the theory above the critical temperature (r > 0) the IR regularization 
is provided by the replacement /c^ — > /c^ + r in the denominator of (A.3)). This becomes 
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impossible for the case r < 0, which we are mostly interested in here. The adequate 
language is then provided by the Legendre transform (effective action) and use of the 
loop expansion or the 1/A^ expansion instead of the primitive perturbation theory. This 
is not convenient, however, for the practical calculation of the renormalization constants. 
Fortunately, in the MS scheme the counterterms are polynomial in IR regulators, and 
the results obtained for them in the region r > can be directly used for r < 0; see 
e.g. the discussion section 3.36 in [21]. Furthermore, the constants Z are independent 
on the specific choice of the IR regularization. From the calculational viewpoints, it is 
more convenient to set r = in the action (and in the propagator (A. 3)) and cut off 
the momentum integrals aX k = m (by dimension, r ~ m'^). Integrals over frequencies 
(or times) are elementary, and one obtains: 

Diagram - 2a /-^ / ^^ /^(A:)D(g)D(|k + q|) 

uiagram - za J ^^^^^ J ^^^^^ ^^^^(k + q)2[A;2 + g2 + (^ + q)2] ' ^^-^^ 

k>m q>m 

where we have denoted 

D{k) = [w{k/iJi)y-' + 1]. (A.6) 

The integrand in (A. 5) depends only on three independent variables: the moduli 
/c = k, g = q and the angle d between the directions k and q, so that kq = kq cosd. 
Thus the expression (A. 5) can be written as a linear combinations of the integrals of 
the form 

dk f dc[ 



(27r)'^ J (27r)'^ 

k>m q>m 



oo 



Sj fdk fdq A:VA:"^g°^|k + q|°3 

-ri+s^uv _L ^^2ri,2 I ^2 I ('V I r.,^21// y^-' ) 



{271^^ J k^+^ J gi+^''(k + q)2[A;2 + g2 + (k + q) 

m m 

where each ai is either 1 or e — y, the brackets mean the averaging over the angle •& 
normalized as ((!)) = 1 and Sd = 2%'^'^ /T(d/2) is the area of the unit sphere in the 
d-dimensional space. 

From dimensional considerations it is obvious that for the integrals I{m) in (A. 7) 
we have 

J(m) = m-2^+"i23 j(^ ^ x)^ ^^23 = Y^ ai. (A.8) 

i 

Thus 



and 



VrJim) = {-2e + ai23)/(m), Vm = mdm (A. 9) 

I{m) = -m-2-+"^^^^--^ -V^I{m)U=^. (A.IO) 

[2e - ai2z) 

Here the pole is isolated explicitly. The expression 'Dml{fn)\m=i is finite at e = y = 0, 
and we can set e = y = in it. Then all these integrals become equal (all ai also become 
0). 
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The factors 772~2e+ai23 ^iii form dimensionless ratios like {m/ ji)^^'^^ or {m/ ji)^^^^ 
with the /i-dependent factors in expressions (A.l) and (A. 6). Since we are interested 
only in the pole parts we will replace such ratios by unities. Thus we have to calculate 

TZ = -Vml{m)\ra=l,s=y=0- (A.H) 

Since m appears in I{m) only in the lower limits of integration, the differentiation gives 

oo 

'^ = J "^^^^ (l + g2 + 2gcos^)(l + g2 + gcos^) ^^- ^^■^'^^ 

1 

Here we have rescaled the variable q/m -^ q, so that new q is dimensionless and m has 

disappeared. The factor 2 (there are two equal contributions since the integrand was 

symmetrical in q and k) cancels with 1/2 from [/c^ + g^ + (k + q)^] — > 2(1 + g^ + g cos i}). 

For e = 0, the prefactor S^/{2'k)'^'^ in (A. 7) is replaced with l/647r^, while the 

angular averaging acquires the form 

2 n 
((...)) = - / sin^t?... . 

TT Jo 

Calculating the resulting double integral gives TZ = (3/2) ln(4/3). 

Now consider the total cofactor which contains the poles in e and y. It comes from 
the denominators in (A. 12) and has the form: 

^ f 1 3^2 3w w^ ] ,, _,^, 

V=< — + - + -, - + -}. A.13 

\2e 2e-2{e-y) 2e - {e - y) 2e - ?>{e - y) \ ^ ' 

Then constant Zi which cancels the poles in (A.l) will have the form 

Zi = l nV (A.14) 

6 

(2(7, being an overall factor in the first term of (A.l) and in (A. 5), does not enter the 

expression for Zf). The anomalous dimension is 

7i = V^ In Zi 

in renormalized variables using the chain rule we obtain {u = A/IGtt^) 

7i = [Pudu + Pwd^oj] In Zi 
with f3 functions from (20). In our approximation, 

V^ = -eV^ + {-y + e)V^ and lnZi = -— 7^P. 

D 

This finally gives 

\2 

71 = — 7^(l + wf = u^{l + w)"^ ln(4/3), (A.15) 

o 

as already stated in (34). For w = one obtains 71 = -u^ ln(4/3), in agreement with the 

result known for the model A; see e.g. [21, 34]. 
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In the same manner we can derive the other results given in (34). For 72,3 the 
integral is quadratic and it should be expanded in the external frequency and momentum 
to Q and p^; the coefficients will be logarithmic integrals of the type (A. 5) and we can 
proceed as before for (A. 7). For 745 this is much simpler because the diagrams are one- 
loop ones, they are logarithmic, there is only one momentum k and the trick involving 
the differentiation P^ is not needed. 
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